Cell organization in soft media due to active mechanosensing 
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Adhering cells actively probe the mechanical properties of their environment and use the result- 
ing information to position and orient themselves. We show that a large body of experimental 
observations can be consistently explained from one unifying principle, namely that cells strengthen 
contacts and cytoskeleton in the direction of large effective stiffness. Using linear elasticity theory 
to model the extracellular environment, we calculate optimal cell organization for several situations 
of interest and find excellent agreement with experiments for fibroblasts, both on elastic substrates 
and in collagen gels: cells orient in the direction of external tensile strain, they orient parallel and 
normal to free and clamped surfaces, respectively, and they interact elastically to form strings. Our 
method can be applied for rational design of tissue equivalents. Moreover our results indicate that 
the concept of contact guidance has to be reevaluated. We also suggest that cell-matrix contacts 
are upregulated by large effective stiffness in the environment because in this way, build-up of force 
is more efficient. 



The mechanical activity of adherent cells usually is at- 
tributed to their physiological function. For example, fi- 
broblasts are believed to maintain the structural integrity 
of connective tissue and to participate in wound healing 
by actively pulling on their environment. During recent 
years, it has become clear that there is another impor- 
tant role for mechanical activity of adherent cells: by 
pulling on their environment, cells can actively sense its 
mechanical properties and react to it in a specific way 
QilMQ- Harris and coworkers observed surprisingly large 
tension fields for fibroblasts on elastic substrates, which 
induce mechanical activity of other cells, even when lo- 
cated at considerable distance 3 . When plated on elastic 
substrates of increased rigidity, many cell types show in- 
creased spreading and better developed stress fibers and 
focal adhesions 5]. Fibroblasts on elastic substrates ori- 
ent in the direction of tensile strain and locomote 
in favor of regions of larger rigidity or tensile strain . 
The same response has been reported for vascular smooth 
muscle cells on rigidity gradients Similar observa- 
tions have been reported numerous times also for tissue 
cells in hydrogels. For fibroblasts in collagen gels. Bell 
and coworkers not only found that traction considerably 
contracts the gel, but also reported orientational effects: 
cells align along the direction of pull between fixed points 
and parallel to free surfaces |^. When a collagen gel is 
stretched uniaxially, cells orient in the direction of prin- 
cipal strain Moreover, cells align in a nose-to-tail 
configuration, thus forming strings running in parallel to 
the direction of external strain. If a collagen gel is cut 
perpendicular to the direction of tensile strain and if cells 
are present in sufficient numbers, they round up and re- 
orient parallel to the free surface introduced llj. 

The response of adherent animal cells to mechanical 
input has evolved in the physiological context of a mul- 
ticellular organism and plays a crucial role in develop- 
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ment, tissue maintenance, angiogenesis, wound contrac- 
tion, inflammation and metastasis. Recently, there has 
been a large experimental effort to understand its molec- 
ular basis. A growing body of evidence suggests that 
focal adhesions based on transmembrane receptors from 
the integrin family act as mechanosensors which directly 
feed into cellular regulation jl2j |. In particular, the ap- 
plication of external force to focal adhesions leads to 
their structural reinforcement and strong signalling activ- 
ity [isl- Q, El and internally generated force correlates 
with the state of aggregation of mature focal adhesions 
IT^ . The exact mechanism of the mechanosensor at 
focal adhesions is still unknown, although structural re- 
organization of the whole aggregate or conformational 
changes of specific molecules are likely candidates. Al- 
though focal adhesions are characteristic for cells cul- 
tured on flat and rigid substrates, cells in a soft envi- 
ronment develop similar cell-matrix contacts which pre- 
sumably have the same mechanosensory function Igj. 

As a result of active mechanosensing at cell-matrix con- 
tacts, cells remodel their contacts and cytoskeleton. In 
particular, they might change position and become ori- 
ented in a certain direction, depending on the mechanical 
properties of their environment. Although cellular be- 
havior in principle results from very complex regulatory 
processes, here we show that the typical cellular reac- 
tion to mechanical input seems to be a simple preference 
for large effective stiffness: starting from this principle, 
we are able to explain many experimental findings which 
have been reported for the behavior of adherent cells both 
on elastic substrates and in hydrogels. In order to make 
these predictions, we have to calculate how stress and 
strain propagates in the extracellular environment. For 
this purpose, we model it with linear elasticity theory 
and solve the elastic equations for different geometries 
and boundary conditions of interest. We then calculate 
the position and orientation in which the cell senses maxi- 
mal effective stiffness in its local environment. Predicting 
cell organization in a soft medium not only contributes 
to a better understanding of many physiological situa- 
tions, but also is of large practical value for application 
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in tissue engineering, e.g. when culturing fibroblasts in 
collagen gels. 

Theory 

Optimization principle for single contact. Mo- 
tivated mainly by recent experiments with elastic sub- 
strates IMS 13 J we suggest that an adherent cell positions 
and orients itself in such a way that it senses maximal ef- 
fective stiffness in its environment. In these experiments, 
the most relevant input for cellular decision making is 
local elasticity of the surrounding environment. Thus 
we first have to calculate how stress and strain in the 
medium is propagated towards the cell. These calcula- 
tions are in general very complicated and will be pre- 
sented below for different situations of interest. In or- 
der to keep our calculations feasible, we assume that the 
extracellular environment is described by isotropic lin- 
ear elasticity theory. This holds true for most synthetic 
elastic substrates and might be a reasonable assump- 
tion for hydrogels. Hence, there are two elastic moduli, 
the Young modulus E (which describes rigidity) and the 
Poisson ratio ly (which describes the relative weight of 
compression and shear modes). In practice, E will be on 
the order of kPa, which is a typical physiological value 
for tissue stiffness. In most situations, v is expected to 
be close to 1/2 (the value for an incompressible medium), 
but other values might be realized in future applications. 

We then ask in which way the cell will organize itself 
if it probes its local environment by actively pulling on 
it. In particular, we aim to define a quantity which de- 
scribes the kind of information which the cell can extract 
from its soft environment with the help of its contractile 
machinery. We first consider a single cell-matrix contact 
and suggest that an appropriate choice is the work W 
which the cell has to invest into the surrounding elastic 
medium in order to build up some force F at the contact 
position fc- As we will show now, the quantity W can be 
used to describe the effects of increased rigidity E and 
prestrain in the elastic environment on an equal basis. 
Therefore ly is a measure for the effective stiffness of the 
elastic environment as probed through a single contact. 

In the absence of prestrain, the work W invested into 
the environment is 

Wo = ^l d^r C,,,iu'r^{r)uli{f) = ^F ■ ^r,), (1) 

where summation over repeated indices is implied. Here 
Uc is the displacement caused by the cell, u'lj (f) the cor- 
responding strain tensor and Cijki the elastic constant 
tensor based on E and i^. The volume integral runs over 
the whole space filled with extracellular material and its 
conversion into a local expression requires partial inte- 
gration and use of the mechanical equilibrium conditions 
(details of our calculations will be published elsewhere). 
Formally, the self-energy of a given contact diverges for a 
point force, but this divergence can easily be removed by 



assuming distributed force. Since displacement decreases 
with increasing rigidity (uij ~ ^/E), the cell has to in- 
vest less work Wq in order to achieve a certain force F 
when rigidity E increases. Hence, the cell senses maximal 
stiffness at the contact when Wq is minimal. 

In a homogeneous medium, the elastic constants do 
not change and Wq is a constant. However, the work 
W needed to build up some force F at the contact po- 
sition fc can vary due to the presence of prestrain. The 
corresponding contribution to W is 

AW = J £t C,jMU%{r)uli{r) = F ■ u^ir,), (2) 

where Ue is the displacement caused by the external 
strain and ufj(r) the corresponding strain tensor. Since a 
negative AW reduces the cellular work W = Wq + AW, 
as does a larger rigidity E, it represents an effective 
stiffening of the environment (strain- stiffening). Corre- 
spondingly, a positive AW represents an effective soften- 
ing with respect to the unstrained medium. Therefore the 
quantity W allows to characterize the local elastic input 
available to an actively mechanosensing cell within the 
unifying concept of effective stiffness, independent of its 
physical origin, which might be rigidity or prestrain. In 
the following, we will identify optimal cell position and 
orientation with the specific force pattern which mini- 
mizes the quantity W . In the sense described here, this 
corresponds to a cellular preference for maximal effective 
stiffness in its local elastic environment. 

It is important to note that conceptually the princi- 
ple suggested here does not imply that the cell actually 
minimizes the work W invested into its soft environment. 
Instead we suggest here that calculating the quantity W 
for different situations of interest is an appropriate mea- 
sure for the kind of information a cell can extract from 
its elastic environment through active mechanosensing. 
The real justification of our model will be its success in 
explaining a large body of experimental data (see results 
section). Nevertheless, below we will also present some 
potential mechanism for the cellular preference for effec- 
tive stiffness, which in fact uses the quantity W not as 
a characterization of the external environment, but as a 
relevant quantity for some internal mechanism. 

Optimization principle for cellular force pat- 
tern. Different contacts are coupled through the actin 
cytoskeleton in such a way that overall force balance is 
ensured. We account for this constraint by considering 
only pairs of opposing forces. In elasticity theory, such 
a pinching force pattern is known as an anisotropic force 
contraction dipole, that is the tensor Pij = PniUj, where 
P is the dipole strength, the product of force magnitude 
and force separation, and n its orientation JJ, 20, 2lj. 
Typical cellular dipoles have been measured to be of the 
order of P « — 10~^^ J (this corresponds to two forces of 
200 nN each, separated by a distance of 60 /xm) [23| . The 
effect of external strain on the work required to build up 
the force dipole Pij at the cell position Tc can be written 
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(a) (b) 

FIG. 1: An adherent cell actively pulls on its soft environment 
through cell-matrix contacts. Experimentally, one finds that 
cells orient themselves in the direction of maximal stiffness of 
the environment. With this cartoon, we present one possible 
mechanism by which active mechanosensing in an elastically 
anisotropic medium might lead to cell orientation. The lo- 
cal elastic environment is represented by linear springs with 
different spring constants K, as indicated by differently sized 
springs. For upregulation of a contact, the cell has to invest 
the work F^/2K. Therefore, upregulation is more efficient 
for larger K. (a) In an isotropic environment, all spring con- 
stants are the same, growth at different contacts is similar and 
the cell does not orient, (b) If spring constants are largest in 
one specific direction, the corresponding contacts outgrow the 
others and the cell orients in the direction of maximal stiff- 
ness of the environment. In this paper, we use the cellular 
preference for large effective stiffness and modelling of the ex- 
tracellular environment by linear elasticity theory to predict 
cell positioning and orientation in soft media. 



as 

AW ^ P.,u%{r,). (3) 

Like in the case of a single contact, optimal cell orga- 
nization can be identified with the specific force pattern 
which minimizes I^W . It follows directly from Eq. (PJ 
that due to the contractile activity of the cell {P < 0), 
tensile strain [u^^ > 0) will always be favorable (negative 
AW). In contrast to compressive strain {ufj < 0), it cor- 
responds to an effective increase in stiffness. Note that 
in contrast to position, orientation and sign, the magni- 
tude |P| of the cellular force pattern does not matter in 
the model presented here. This reflects the fact that here 
we aim to characterize the mechanical properties of the 
extracellular environment sensed by the cell, rather than 
the process of active mechanosensing itself. 

Possible origin of optimization principle. Our 
modelling starts from the phenomenological observation 
that cells seem to prefer maximal effective stiffness in 
their environment. Although it can be justified by its 
large success in explaining experimental observations (see 
results section), we also want to suggest a possible mech- 
anism for our main assumption. For this purpose, we use 
a simple one-dimensional analogue. Consider the extra- 
cellular environment to act like a linear spring with spring 
constant K, on which the cell is pulling through a single 
cell-matrix contact. Recent experiments on focal adhe- 
sions ^3 suggest that upregulation of contact growth 



is related to reaching a certain threshold in force F, al- 
though the details of how force affects regulation are just 
beginning to emerge 12]. In order to build up sufficiently 
large force F, the cell has to invest energy W — F'^/2K 
into the spring. Thus, the stiffer the spring (the larger 
K), the less work is needed and the more efficient the 
build-up of force will be. An equivalent viewpoint is to 
assume that the cell invests the power L into stretching 
the spring. Then it takes the time t = F^ /2KL to reach 
the force F. Therefore, a specific contact will grow faster 
than the other contacts if it encounters a larger stiffness 
K. In principle, the cellular program could also be geared 
towards achieving a certain displacement of the surround- 
ing material, which would result in a preference for effec- 
tive softness of the environment. However, this scenario 
would imply the existence of some additional mechanism 
for outside-in signaling. It is more realistic to consider 
that force activates the cellular response through a cer- 
tain displacement of elastic components located inside 
the cell. Since internal displacement and force are ex- 
pected to be linearly related through another (internal) 
spring constant, one arrives at the same result. 

Adhering cells probe the mechanical properties of their 
environment by pulling at many cell-matrix contacts si- 
multaneously (Fig.P). At each newly formed contact the 
cell is expected to pull with a similar investment of re- 
sources (e.g. actin, myosin or ATP). However, each con- 
tact encounters a different elastic environment, each of 
which can be represented by a different spring constant. 
In an isotropic situation (Fig. ^), all spring constants 
are equal, the contacts have similar growth behavior and 
there is no reason for a cell to orient. Experimentally, the 
cell adopts a round or stellate morphology, depending on 
the number of cell-matrix contacts. In an anisotropic 
situation (Fig. ^), build-up of force is more efficient in 
one specific direction and contacts in this direction will 
eventually outgrow the other ones. Here, the anisotropic 
elastic properties of the medium provide an orientational 
clue for the cell, which orients along the direction of max- 
imal effective stiffness. Depending on e.g. the presence of 
motility factors, this orientation response might be fol- 
lowed by cell locomotion. 

Results 

Homogeneous external strain. We first consider a 
cell interacting with homogeneous external strain, either 
on the top surface of a rectangular slab of elastic ma- 
terial (elastic substrate) or inside an infinite elastic ma- 
terial (hydrogel). In both cases, the equations of three- 
dimensional isotropic elasticity give 

AW ^ -—[{l + iy)cos'^e -v] (4) 
E 

where 9 is the orientation angle relative to the direction 
of the externally applied tensile stress p < 0. Optimal 
cell orientation corresponds to minimal AW, which is 
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FIG. 2: Adjusting cell position and orientation in such a way 
that the cell senses maximal effective stiffness in its environ- 
ment is equivalent to minimizing the quantity W, the amount 
of work the cell invests into the elastic surroundings in the 
presence of external strain. In the presence of mechanical 
activity, sample boundaries induce external strain which can 
result in different cell organization, (a) AW for a cell with 
dipole strength P which is a distance d away from the sur- 
face of an elastic halfspace with rigidity E, plotted in units 
of P^/Ed^ as a function of angle 6 between cell orientation 
and surface normal (rescaled by 2567r). Solid and dashed lines 
correspond to Poisson ratios i/ = 1/2 and v = 0, respectively. 
Irrespective of v, the optimal orientations (minimal AW) are 
perpendicular {9 = 0) and parallel {6 = 7r/2) to the surface 
for clamped and free boundaries, respectively. Since |Aiy| 
increases if d decreases, the overall mechanical activity of a 
cell increases towards a clamped surface {AW < 0), but de- 
creases towards a free surface (AW > 0). (b) AW for a cell 
in an elastic sphere of radius R, plotted in units of P/ER^ 
as a function of distance r to the sphere center in units of R 
for f = 1/3 (rescaled by 15/8). Solid and dashed lines are 
parallel {9 = 7r/2) and perpendicular {9 — 0) orientations, 
respectively (all other orientations yields curves which lie in- 
between the ones shown). Like in an elastic halfspace, paral- 
lel and perpendicular orientations are favored (minimal AW) 
for free and clamped boundaries, respectively. For clamped 
boundaries, mechanical activity is favored (smaller AW) to- 
wards the surface. For free boundaries, mechanical activity is 
disfavored (larger AW) towards the surface. 



achieved for 9 — 0, irrespective of the Poisson ratio v. 
Thus the cell orients preferentially with the direction of 
stretch. This behavior is indeed observed experimentally, 
both for fibroblast on elastic substrates Q and in colla- 
gen gels 0, 0. Since AW decreases with increasing 
rigidity E, the elastic effects discussed here will be ob- 
served only in a soft environment, namely with rigidity E 
around kPa, which is a typical physiological value for tis- 
sue stiffness. For stiffer substrates the variations in AW 
for different contact positions might become too small to 
induce an orientation response. 

Boundaries. In a physiological context, cells are of- 
ten close to boundaries, like the surface of a tissue or 
organ. In the presence of cell traction, boundaries alter 
the strain with respect to a homogeneous infinite medium 
by a boundary-induced strain {image strain), which has 
a similar effect like external strain. In this way, cells can 
actively sense not only the presence of a close-by surface. 




FIG. 3: Predicted cell orientation in a hydrogel close to a 
surface (a,b) and on elastic substrates (c,d). (a) Cells prefer 
the direction of maximal effective stiffness. Thus, they orient 
perpendicular to a clamped surface, (b) For a free surface, 
this direction is parallel to the surface, (c) Cells close to a 
boundary between soft (left) and rigid (right) regions prefer 
analogous orientations as cells close to clamped and free sur- 
faces in a hydrogel, respectively, (d) Cells interact elastically 
to form strings, because in nose-to-tail alignment, the me- 
chanical activity of one cell triggers the one of the other cell, 
thereby forming a positive feedback loop. 



but also its shape and boundary conditions. In order 
to predict the effect of boundaries on cell organization, 
we study a semi-infinite space with a planar surface, for 
which the elastic equations can be solved exactly |23l |. 
The details of the boundary conditions in a physiological 
context can be very complicated. Here we address two 
fundamental reference cases, namely free and clamped 
boundary conditions, for which normal stress and dis- 
placement, respectively, vanish at the surface. Consider 
a force dipole which is a distance d away from the planar 
surface and has an angle of orientation 9 relative to the 
surface normal. We find 

p2 

AW ^ —{a^ + Kcos'^9 + c^cos'^9), (5) 
Ea-^ 

where the coefficients a^, , bu and Cu are complicated func- 
tions of v. AW scales quadratically with P because 
the image strain scales linearly with P (in other words, 
the force dipole interacts with its images). For free and 
clamped surfaces, the coefficients are positive and nega- 
tive, respectively, irrespective of v. Therefore, the opti- 
mal configurations (minimal AW) are parallel {9 = 7r/2) 
and perpendicular {9 — 0) for free and clamped bound- 
aries, respectively, as plotted in Fig. Et- A schematic 
representation (Fig. 13^ and b) provides a simple inter- 
pretation: for clamped (free) boundary conditions, the 
cell senses maximal stiffness perpendicular (parallel) to 
the surface. One may think of a clamped (free) surface 
as the interface between the medium and a imaginary 
medium of infinite (vanishing) rigidity, which effectively 
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rigidifies (softens) the medium towards the boundary. In 
general, we find that free and clamped boundary condi- 
tions always have opposite efi'ects, albeit with one essen- 
tial difi^erence: for clamped boundaries, mechanical ac- 
tivity of cells is favored and cells can amplify this effect 
by adjusting orientation. For free boundaries, mechanical 
activity of cells is disfavored and the orientation response 
is an aversion response. 

Experimentally, it is well known that mechanical ac- 
tivity of cells increases for clamped boundary conditions 
[23 | . The predicted orientation effects close to boundaries 
have been observed numerous times, e.g. the parallel ori- 
entation of cells close to free surfaces 0. Our model 
predicts the same orientation effects for an elastic sub- 
strate with two regions of different rigidities (Fig. 
cells on the soft and stiff sides of the boundary orient 
perpendicular and parallel to it, respectively. Indeed fi- 
broblasts migrating from a soft to a stiff region keep their 
perpendicular orientation and cross over to the stiff side, 
while fibroblasts migrating from a stiff to a soft region 
do not cross the boundary, but turn by 90 degrees and 
move parallel to the boundary 0. 

Finite sized sample. As an example for a finite sized 
sample, we consider the elastic sphere with radius R. The 
elastic equations can be solved exactly by using an ex- 
pansion in terms of vector spherical harmonics |2a|. We 
find 
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where r denotes distance to the sphere center, 9 the ori- 
entation in respect to the radial direction and fi, is an 
infinite sum over all angular momenta, which does not 
change qualitatively as v is varied. In regard to orienta- 
tion, we find the same results as for the elastic halfspace 
(compare Fig. : cells will orient parallel (perpendicu- 
lar) to free (clamped) surfaces, respectively. We also find 
a similar result for the effect of distance to the surface: 
for free (clamped) boundary conditions, a small (large) 
distance to the sphere center is more favorable, since the 
surface favors (disfavors) mechanical activity. The new 
aspect here is the role of sphere radius R: since \AW\ 
increases when sphere radius R decreases, one can ef- 
fectively rigidify (soften) a material with clamped (free) 
boundaries by reducing system size. Our predictions 
could be tested using e.g. fibroblast-populated collagen 
microspheres, an assay which has been introduced to 
stud y c ompaction of tissue equivalents at high cell den- 
sity [26] . Since here we are mainly concerned with single 
cell effects, we suggest to modify this assay in such a way 
as to monitor the organization of isolated cells close to 
the sphere surface at low cell density and as a function 
of varying sphere radius. 

Cooperative effects. Up to now we have been dis- 
cussing single cell effects and we now turn to coopera- 
tive effects. In particular, we now consider the case that 
external strain is caused by the traction of other cells, 
which amounts to an elastic interaction of cells. Even if 




FIG. 4: Monte Carlo simulations of elastically interacting cells 
in an external strain field. The temperature used in the sim- 
ulation represents the stochastic element of the process of cell 
organization. Without external strain, cells form strings. In 
its presence, strings align in parallel. 



all cells initially have isotropic force patterns, they will 
sense anisotropic strain fields and start to orient them- 
selves. For the simplest case of two cells, we find 



9.(01,02,0), 



(7) 



where r is the distance between the force dipoles and g^, 
is a complicated function of v and the three orientational 
degrees of freedom. Our calculation shows that AW has 
a pronounced minimum for completely aligned dipoles, 
independent of (Fig. |3i). In this configuration, both 
cells sense maximal effective stiffness, because maximal 
strain stiffening occurs along the axis of contraction. This 
finding suggests that a common pattern for the organiza- 
tion of elastically interacting cells will be the formation 
of strings of cells. They might close into rings, so that 
each cell can be fully activated by its two neighbors. We 
used Monte Carlo simulations to obtain a typical configu- 
ration of elastically interacting cells in an external strain 
field (Fig. 01). The temperature of the Monte Carlo simu- 
lation represents the stochastic nature of the orientation 
response. We find strings of cells aligned in parallel with 
the external strain, exactly as observed experimentally 

IS-.. 

It is important to note that there is a positive feedback 
for cell alignment: the more cells orient in one direction, 
the stronger becomes the input for other cells to adopt 
the same orientation. For example, it has been reported 
that when boundary condition are changed from clamped 
to free by cutting the collagen gel, fibroblasts show the 
predicted reorientation by 90 deg rees only when suffi- 
ciently many cells are present In fact our calcula- 
tions and simulations show that cells can orient in par- 
allel even with respect to clamped boundaries, if there 
are sufficiently many cells such that the direct elastic in- 
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teraction between cells dominates the single cell response 
of perpendicular orientation. In practise, the single cell 
response might also be disturbed because elastic signals 
could be screened by traction of randomly oriented cells. 
Indeed such an effect has been reported for experiments 
with elastic substrates Q. 



Discussion 

It has long been known, especially in the medical and 
bioengineering communities, that cell organization in soft 
media is strongly influenced by the mechanical properties 
of the environment. Here we presented a model which is 
able to explain numerous experimental observations that 
have been reported for organization of cells (especially 
fibroblasts) both on elastic substrates and in hydrogels. 
The excellent agreement of our results with experiments 
suggests that cell organization can be predicted from lo- 
cal mechanical properties which the cell actively senses 
in its environment. In fact the only property of cellular 
regulation which enters our model is the assumption that 
cells locally prefer large effective stiffness. Otherwise our 
modelling focuses on the elastic properties of the extra- 
cellular environment. 

Modelling the soft environment of cells as an isotropic 
elastic medium is certainly a good assumption for elas- 
tic substrates. The situation is more complicated for 
hydrogels, in particular because they might not behave 
elastically and because they feature fiber degrees of free- 
dom. Cell organization in gels is often explained by con- 
tact guidance, the alignment of cells along topographic 
features like collagen fibers. Since fibers can become 
aligned due to cell traction, contact guidance provides 
a long-ranged and persistent mechanism for cellular self- 
organization in tissue equivalents This process has 
been modeled before. In the theory of Ref. j23|, flux 
equations for cellular and matrix densities are combined 
with mechanical equations which include cells as centers 
of isotropic contraction. This might be a good model 
for chondrocytes, which tend to show a spherical mor- 
phology. The anisotropic biphasic theory (ABT) from 
Ref. [23 aims at cells like fibroblasts and smooth muscle 
cells, whose typical morphology in tissue equivalents is 
bipolar. ABT introduces a cell orientation tensor, which 
is coupled to a fiber orientation tensor, since cells are as- 
sumed to react foremost to fiber degrees of freedom. In 
our model, the force dipole tensor represents cell orien- 
tation as does the cell orientation tensor in ABT, but it 
is coupled to elastic degrees of freedom, since cells are 
assumed to react foremost to large effective stiffness. 

Because models for contact guidance in tissue equiv- 
alents focus on fiber degrees of freedom and high cell 
densities, they do not explain the single cell responses 
observed on elastic substrates, where contact guidance 
usually is ruled out Q- The large predictive power 
of our model for elastic substrate experiments suggests 
that active mechanosensing by single cells might also be 



involved with cell organization in hydrogels. However, 
for the collagen assay from Ref. (2^J it has been shown 
that as a result of external strain, fibers become rear- 
ranged and stress relaxes towards zero. In a matrix which 
cannot support any stress, our elastic considerations do 
not apply and contact guidance through formerly aligned 
fibers might be the only relevant clue for cell organiza- 
tion pQj . However, it is important to note that in our 
model, stress is actively generated by cells and thus needs 
to be supported only over time scales in which the cell 
actively senses the mechanical properties of its environ- 
ment. In particular, if fiber alignment has resulted in 
some anisotropic mechanical environment, the cell might 
sense the anisotropic mechanical properties of the matrix 
and orient itself correspondingly. This might explain why 
cells have been found to align to a greater extent with 
respect to external strain than the surrounding collagen 
fibrils '2^ and why our modelling is also successful for 
hydrogels. In general, future experiments are needed to 
clarify the relative importance of topographic versus me- 
chanical clues for cell organization in hydrogels, while 
future modelling is needed to account for the mechanical 
(in particular, viscoelastic) properties of hydrogels. 

We also want to point out that contact guidance is 
a bidirectional clue and provides only guidance, in con- 
trast to external elastic strain, which provides taxis. In 
our model, taxis is reflected by the position dependence 
of AW. For example, our theory not only predicts that 
cells prefer to orient parallel to free boundaries, but also 
that cells prefer to move away from them. Moreover a 
simple preference for cell alignment along fibers does not 
predict what cells do if they encounter a fiber junction 
in the gel. Although we are not concerned with cell loco- 
motion here, our modelling would suggest that cells pre- 
fer the fiber under largest tension, exactly as has been 
observed experimentally for neutrophils migrating in hu- 
man amnion |30| . 

During recent years, the regulated response to mechan- 
ical input by single cells has been studied experimen- 
tally in larger detail. There is a growing body of ev- 
idence now that integrin-based cell-matrix contacts act 
as local mechanosensors which channel mechanical in- 
formation about the environment directly into cellular 
decision making. Although this does not concern our 
modelling directly, here we suggested that the upregula- 
tion of growth of cell-matrix contacts in a stiff environ- 
ment might result from the fact that it is triggered by 
a threshold in force, whose build-up is more efficient for 
larger stiffness. An equivalent viewpoint is that growth 
of cell-matrix contacts is faster on stiffer substrates. As 
experimental test for this hypothesis, we suggest corre- 
lation studies for growth of cell-matrix contacts and cel- 
lular organization, especially close to sample boundaries, 
where cells can amplify the mechanical input provided by 
boundary induced strain through active mechanosensing. 
Quantitative data about growth behavior of cell matrix 
contact will allow us to further refine our model in a more 
quantitative way, possibly also including modelling of eel- 
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lular features like morphology and force pattern, which 
are not the focus of this work. 

It is important to note that our model suggests com- 
pletely different behavior for cells then one would expect 
for physically inert particles interacting with a soft ma- 
trix (like mobile inclusions in a metal). Although the 
expression in Eq. (j2Jl is similar to the interaction poten- 
tial for the physical system of force dipoles in an external 
strain field, it has the opposite sign, because in the phys- 
ical system, one has to minimize the composite energy 
of defect and medium 19, 20] . The physical potential 
has been used before to model elastic interactions of cells 
without any regulatory response and has been shown to 
lead to aggregation behavior similar to the one of elec- 
tric quadrupoles (2ll |. while the model introduced here 
leads to aggregation behavior similar to the one of elec- 
tric dipoles |3ll| . 

Since elastic effects are long-ranged and propagate 
quickly, they provide an appealing mechanism for signal 
transduction for mechanically active cells in soft media. 
However, they are also unspecific and cells might not 
be able to distinguish between different sources. On the 
other hand, additional information channels, like soluble 
ligands, will certainly supplement elastic signals. More- 
over, cells in highly differentiated organisms are likely to 



interprete mechanical signals only in their own physio- 
logical context, which is more restricted than for cells in 
an arbitrary environment. 

In summary, we have presented an optimization prin- 
ciple in linear elasticity theory that allows to predict cell 
organization in soft media in excellent agreement with 
a large body of experiments. Moreover we have sug- 
gested a mechanism which links the cellular preference 
for large effective stiffness to growth of cell-matrix con- 
tacts. Our modelling results in many interesting predic- 
tions which now can be checked experimentally. In the 
future, our model might be extended to high cell densi- 
ties and strong cooperative effects, which are character- 
istic for tissue equivalents. We expect that then it can 
be used for rational design in tissue engineering. For ex- 
ample, using numerical (finite element) methods, one can 
use it to optimize protocols for the design of tissue equiv- 
alents for implants in regard to geometry and boundary 
conditions. 
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